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Abstract 

We study cosmic-ray anomaly observed by PAMELA based on Ee inspired extra U(l) model with 
54 flavor symmetry. In our model, the lightest flavon has very long lifetime of O(10 18 ) second which is 
longer than the age of the universe, but not long enough to explain the PAMELA result 0(1O 26 ) sec. 
Such a situation could be avoidable by considering that the flavon is not the dominant component of 
dark matters. However non-thermalizing the flavon is needed to obtain proper relic density. This relates 
reheating temperature of the universe with seesaw mass scale. If we assume this flavon is a particle 
decaying into positron (or electron), the seesaw mass scale is constrained by reheating temperature. 
Thus we find an interesting result that the allowed region is around O(10 12 ) GeV, which is consistent 
with our original result. 
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1 Introduction 



It is one of the important task to build more economical models in (non-)Abelian flavor symmetries. In such 
the framework of models with non-renormalizable operators especially, couplings of the terms are usually 
suppressed by high energy cut-off scale. Therefore gauge singlet neutral bosons which couple to the term 
may play an important role not only to construct the mass matrix forms but also to be promising dark 
matter candidates, because they could be less-interactive enough. It is recently known that the dark matter 
can be a good candidate to explain PAMELA data [1]. Subsequently, there are many attempts to explain the 
positron anomaly by annihilation [2] or decay of the dark matter [3] . According to constraint from diffuse 
gamma ray [4], an interpretation by annihilation is almost excluded. Thus the PAMELA result is in favor 
of the decaying dark matter, when it has the life time of F -1 <~ 0(1O 26 ) sec. This is much longer than the 
age of the universe. 

In this paper, we study such a cosmic-ray excess by a singlet scalar (flavon) in S4 flavor model 1 of a 
supersymmetric [5] with E e inspired extra U(l) gauge symmetry [6]. The flavor symmetry is broken by 
vacuum expectation value (VEV) of the flavon and this VEV gives large Majorana masses for right handed 
neutrinos (RHNs). As the flavon couples to standard model particles only through non-renormalizable 
operators, the life time of the flavon could be longer than the age of the universe. However the life time 
of the flavon in our model is not longer than O(I0 26 ) sec. Therefore we consider that the flavon is not 
the dominant component of dark matters. Then the most leading interaction which the flavon has is extra 
U(l)2 interaction and the interaction is extremely weak because it is suppressed by the large mass scale 
of the U(l)z g au ge boson. As the result, annihilation cross section is too small to obtain proper relic 
density of the flavon as long as we assume that the flavon is in thermal equilibrium. So we consider that 
the flavon is never in thermal equilibrium. Non-thermalizing the flavon constrains reheating temperature 
and also relates to right-handed neutrino mass scale. If we assume reheating temperature is constrained 
as 10 7 GeV > Trh > 10 4 GcV, right-handed neutrino mass scale should be around I0 12 GeV, which is 
consistent with our original result. 

The paper is organized as follows. In section 2, we review the basic structure of S4 flavor symmetric 
extra U(l) model. We evaluate density parameter of lightest flavon from PAMELA constraint in section 3, 
and estimate required right-handed neutrino mass scale for explaining it in section 4. Finally we make a 
brief summary in section 5. 



2 The Extra U(l) Model with S4 Flavor Symmetry 
2.1 The Extra U(l) Model 

The basic structure of the extra U(l) model is given as follows [7]. At high energy scale, the gauge symmetry 
of model has two extra U(l)s, which consists maximal subgroup of E 6 as G 2 — Gsm x U(1)x x U(1)z C E e . 
MSSM superfields and additional superfields are embedded in three 27 multiplets of E 6 to cancel anomalies, 
as 27 D {Q, U c , E c , D c , L, N c , H D ,g c , H u , g, S}, where N c are right-handed neutrinos (RHN), g and g c arc 
exotic quarks (g-quark), and S are Gsm singlets, which is illustrated in Table 1. We introduce Gsm x U(1)x 
singlets $ and $ c which develop the intermediate scale VEVs along the D-flat direction of ($) = (& c ), then 
the U(l)z is broken and the RHNs obtain the mass terms. After the symmetry is broken, as the R-parity 
symmetry remains unbroken, G\ = Gsm x U(1)x x R survives at low energy. This is the symmetry of the 
low energy extra U(l) model. 

Within the renormalizable operators, G2 symmetric superpotential is given as follows: 

W 2 = Wo + Ws + Wb, (1) 

W = Y U H U QU C + Y D H D QD C + Y E H D LE C + Y N H U LN C + Y M <$>N C N C , (2) 

W s = kSgg c + XSH U H D , (3) 

W B = X 1 QQg + X 2 g c U c D c + X 3 gE c U c + X 4 g c LQ + X 5 gD c N c . (4) 

Where Wq is the same as the superpotential of the MSSM with the RHNs besides the absence of /i-term, 
and Ws and Wb are the new interactions. In Ws, kSgg c drives the soft SUSY breaking scalar squared 
mass of S to negative through the renormalization group equations (RGEs) and then breaks U(l)x and 

1 It requires exotic scalar quarks, which induce proton decay. In our original work [7], we have shown that proton decay is 
suppressed by 54 flavor symmetry very well. 
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generates mass terms of g-quarks, and XSH U H D is source of the effective /x-term. Therefore, W and W$ 
are phenomenologically necessary. In contrast, Wb leads to very rapid proton decay and must be forbidden. 
This is done by S4 flavor symmetry. 
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Table 1: G2 assignment of fields. Where the x, y and z are charges of U(l)x, U(l)y and U(l)z, and Y is 
hypercharge. 



2.2 S4 Flavor Symmetry 

Non-Abelian group S4 has two singlet representations 1, 1', one doublet representation 2 and two triplet 
representations 3, 3', where 1 is the trivial representation [8]. The most essential structure of £4 group is 
that multiplication of two doublets does not contain triplets. With this property, if g and g c are assigned to 
triplets and the others are assigned to singlets or doublets, then Wb is forbidden. 

The absence of Wb makes g-quarks and proton stable, but the existence of g-quarks which have life time 
longer than 0.1 second spoils the success of Big Ban nucleosynthesis. In order to evade this problem, we 
assign $ c as triplet of S4 and add the non-renormalizable terms: 



w NRB = i$$ c (QQg + g c u c D c + g E c u c + 9 c lq + 9 d c n c ) . 



When $ c develops VEV with 



10" 



(5) 



(6) 



the phcnomenological constraints on the life times of proton and g-quarks are satisfied at the same time [9] , 
and the right-handed neutrino mass scale can be predicted as Mr ~ ($) ~ 10 _6 Mp ~ 10 12 GeV. 

As the flavons <I> and <I> C which are the triggers of flavor violation do not have renormalizable interactions 
with light particles, the lightest flavon [CF) has very long life time. In following sections, we consider 
whether this particle explains PAMELA observation. 



3 Flavon Decay Width 

With the assignment that $ c is S4 triplet and $ is singlet or doublet, the leading term of flavon supcrpotential 
is given by 



W® = — $ 2 ($ c ) 2 
M P 



Solving the potential minimum conditions, we get 

V = |$| = |$ c | ~ (msusYM P /a)i ~ 10 n O -i (^f^) ' GcV, 



(7) 



(8) 



where a <~ 10 2 is required from Eq.(6). Integrating out the heavy RHNs, we get effective seesaw operators 
as follows 



W, 



eff 



1 

Y M <S> 



{Y N H u Lf. 



(9) 
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Here we redefine flavon as perturbation around VEV as $ — » V + $, then Eq.(9) is rewritten as 



Weff = y^y{Y N H^Lf - _*_(r^L)». (10) 



If sleptons, squarks, g-quarks and scalar g-quarks are heavier than CT, this operator gives dominant contri- 
bution to decay width of CF through 



where we assume 



£e// = (11) 



m^^l, (H°)=* (12) 



From the interaction, monochromatic intense diffuse neutrino flux is expected in cosmic-ray at E v = m$/2. 
It is an important signature of the decaying dark matter model. Considering the superpotential W$, it is 
easily shown that only one linear combination of six flavons $1,2,3, $1 2.3 nas super heavy mass around V and 
another five flavons have O(Msusy) masses. As is pointed in ref. [10], this interaction is good candidate 
for explaining PAMELA phenomena, because if V is around 10 16 GeV, then the partial decay width of 
$ — > vH~e + is given by 

r-i(* - vH-e+) = ^|) _1 = 3-5 x 10- sec, (13) 

where 

m^O.leV, v — 246 GeV, = 3 TcV, V = 10 16 GeV. (14) 

The result of Eq.(13) is in good agreement with ref. [3]. However, as V <~ 10 12 GeV in our model, the life 
time of CT is not 10 26 sec but 10 18 sec. So we assume this CF is not the dominant component of dark 
matter (Hip <C SIbm)- Introducing mixing parameter e defined as 

$ = e$ LF + ---, (15) 

where &lf is the lightest flavon field and the dots • • • means contributions from heavier flavons, the partial 
decay width of $lf — > vH~e + is given by 



r _1 ($LF -> ^i?"e + ) = 3.5 x 10 18 e" 2 ^- 



V 2 



sec. (16) 



,10 12 GeV 

In order to explain positron flux observed by PAMELA, density parameter of CF should be 

^ = r ~ 1( **Lr g ~ e+ W = 3.9 x 10-e- 2 f -^^j 2 M (17) 



10 26 sec ' VlO^GeV, 

where we use the observed value of density parameter of dark matter in [11] as 

HflM = O.II//1 2 . (18) 

Before calculating density parameter of CF in next section, we give some comments about possible decay 
channels of flavons. The interactions between flavons are described by 



(m p ) 



2 

3 



V d $i$j$ fc - (10~ 4 GeV)$i*j$ fc , (19) 



from which decay width of heavier flavon to lighter flavons is given by 

F($, -> $j$ fc ) - 10 12 sec" 1 . (20) 
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With this interactions, all heavier flavons decay into LT finally. If LT is heavier than sneutrinos, the 
following operators coming from Eq.(lO) may contribute to CF decay: 

£*nn ~ msusY-^-^LFNN. (21) 

Where N is sneutrino. If CF is heavier than scalar g-quarks, the following operators coming from Eq.(5) 
may contribute too: 

V 

M2* LF 9M> ( 22 ) 

where g is scalar g-quark and q is doublet quark. The contributions from these operators are the same 
order as the contribution from Eq.(ll). Hereafter, we assume CF is lighter than those scalar particles and 
g-quarks. In this case, CT decay channels to these particles are kinematically closed. 



4 Estimation of RHN Mass Scale 

At high temperature T 3> msusY, dominant contribution to \X ~^ come from U(l)z gauge interaction, 
where \ is the particle in thermal bath. This interaction is in thermal equilibrium when 

T( XX -+ $$) ~ ^ > H ~ — (23) 
V Trip 

is satisfied. Note that M P = 2.43 x 10 18 GcV is reduced Planck mass and m P = 1.22 x 10 19 GeV is Planck 
mass. Evading overproduction of gravitino, T < 10 7 GeV must be satisfied with [12]. Then the constraint 
Eq.(23) leads following condition: 

10 10 GcV > V. (24) 

However it is difficult to satisfy this condition in our model. Moreover, if flavon is in thermal equilib- 
rium, there is a problem of over-production of flavon because U(l)z gauge interaction is too weak at low 
temperature to cause appropriate pair annihilation of flavons. Based on this discussion, we assume U(l)z 
gauge interaction is never in thermal equilibrium and the universe starts with low flavon number density 
n$(T RH ) = 0. 

The interaction between U(l)z gauge boson and chiral multiplet (if), ^F), where ip is fermion and \1/ is 
boson, is given by 

Cgauge = iQzA^ ^ Z i [^i.L^i.L + *i^*J ~ * Jfytf i , (25) 

i 

from which thermally averaged cross sections are given by 

X>W-^M> n $ - 55.1CT 8 , (26) 

Y,^^M) n % = 96.1CT 8 , (27) 

$^W->**M>n$ = 96.1CT 8 , (28) 

■4, 

^(a^^^vDnl = 248. 8CT 8 , (29) 



(2tt) 5 V M| 

M g = 2g z zzV = l6g z V, (31) 

where ^ does not contain flavon $ and its fermion partner <f> because n$ and are negligible. Here we 
define the CT production rate Nlf which means how many CFs, are generated per one degree of freedom 
of flavon and its fermion partner. Because the total degrees of freedom is 20, Nlf is bounded from below 



4 



by 1/20 2 . If we require that the masses of flavons and those fermion partners are smaller than 9 TeV, N LF 
is bounded from above by 9 TeV/m F F = 3. Therefore allowed range of Nlf is given by 



3 > N LF > 0.05. 



(32) 



The model dependence of Savon sector affects this analysis only through e and N F p- Using Nlf, the 
Boltzmann equation for the LT number density Ulf- 

n LF = 5A r LF (n^ + n$), (33) 

is given by 

n LF + 3Hn LF = 2480N LF CT 8 , (34) 
where Hubble constant H is given by 

H = l.mJgl , (35) 

mp 

.9, = 341.25. (36) 
Solving Eq.(34) with boundary condition ulf{Trh) = 0, CF to entropy ratio is calculated as follow: 

(n LF \ _ 15 x 2480 x 21m P N LF ( g 2 ^ \ 2 3 

I a Jo 2ttV x 30.67(2tt) 5 \ ) RH ' [ ' 

Note that Eq.(37) does not depend on the definition of U(l)z gauge coupling constant g z and charge 
normalization. Finally wc require the density parameter 



LF 



m LF n L F 

Pc 



satisfies Eq.(17) then we get 



where 



V 



10 12 GcV 



(e 2 N LF ) 



Trh 
10 5 GeV 



s = 2890/cm 3 , 

p c = 1.05 x I0 4 h 2 eV/cm 3 , 



are used [11]. 
If wc assume 



10 7 GcV > T RH > 10 4 GcV > m LF , 
and the model dependent parameters as Eq.(32) and 

1 > e > 0.1, 

we get the allowed range of V as follow: 

V 



12 > 



10 12 GcV 



> 0.09. 



(38) 



(39) 



(40) 
(41) 



(42) 
(43) 

(44) 



This prediction is consistent with our previous result V ~ 10 12 GeV which comes from phenomenological 
constraints for proton life time and g-quark life time [7]. If we use V = 10 12 GeV, reheating temperature is 
predicted as 



T RH = (0.7 - 5.8) x 10 5 GeV. 



(45) 



2 Note that 4 of 24 = 4 X 6 are super-heavy or eaten by gauge boson, where 6 is total number of fiavon supcrficlds (<J>, <I> C ). 
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Finally we give some comments. In the case that CT is dominant component of dark matter as fl^p ~ 
SlflM, which is realized by putting e(10 12 GeV/V) ~2x 10~ 4 , e is bounded from above by 

e < 1.2 x 1CT 4 . (46) 

However, it seems difficult to explain such a small mixing angle. 

If the final states of CF decay contain quarks, anti-proton flux may be significantly induced, which 
conflicts with cosmic-ray observations [13]. This problem is avoided when LT decay mainly into charged 
Higgs which does not couple to quarks. Such a type of model [14] can be constructed with S4 flavor symmetry. 
For example, if all quarks are assigned to 54-singlets and SU{2)w Higgs doublets H u and H D are assigned 
to one S*4-doublet and one S^-singlet respectively, Yukawa interactions between ^-doublet Higgs and quarks 
are forbidden. This assignment also solves Higgs-FCNC problem, as pointed in ref. [15]. However, as there 
is no room for model building in weak interaction, we can not suppress anti-proton production by weak 
interaction such as $ — > vW~e + — > p. It is not obvious whether the effects of weak interaction spoil this 
solution or not, which is left for future work. 

5 Summary 

In this paper, we have considered cosmic-ray anomaly observed by PAMELA based on 5*4 flavor symmetric 
extra U(l) model. Identifying the particle decaying into positron as the lightest flavon which is the sub- 
dominant component of dark matter, reheating temperature of the universe and the mass scale of right 
handed neutrinos were related, and we estimated the mass scale of right handed neutrinos from the relation. 
As a result, the constraint that reheating temperature is low enough to suppress gravitino production gives 
vacuum expectation value of flavon around 10 12 GeV, which is consistent with the prediction of our model. 
This result supports the idea of neutrino seesaw mechanism based on heavy right handed neutrino with 
O(10 12 ) GeV mass. Finally, monochromatic diffuse neutrino flux is expected at E„ = m$/2 in cosmic-ray 
as a predction of the model. 
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